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ABSTRACT 

In  this  paper  we  propose  a  fast  reconstruction  method 
for  a  progressive  subband-decomposed  signal  coding  sys¬ 
tem,.  It  is  shown  that  unlike  the  normal  approach  which 
contains  a  fixed  computational  complexity,  the  compu¬ 
tational  complexity  of  the  proposed  approach  is  propor¬ 
tional  to  the  number  of  refined  coefficients.  Therefore, 
using  the  proposed  approach  in  image  coding  applica¬ 
tions,  we  can  update  the  image  after  receiving  each  new 
coefficient  and  create  a  continuously  refined  perception. 
This  can  be  done  without  any  extra  computational  cost 
compared  to  the  normal  case  where  the  image  is  recon¬ 
structed  after  receiving  a  predefined  number  of  bits. 

1. INTRODUCTION 

Progressive  transmission  of  signals  is  a  mechanism  by 
which  the  encoder’s  output  is  transmitted  in  groups  of 
bits  (packets)  and  the  decoder  produces  a  higher  qual¬ 
ity  replica  of  the  signal  based  on  receiving  each  new 
packet.  Progressive  transmission  is  a  desirable  feature 
in  many  practical  signal  transmission  situations  such  as 
telebrowsing  and  database  retrieval.  Progressive  trans¬ 
mission  also  provides  the  opportunity  for  interrupting 
transmission  when  the  quality  of  the  received  signal  has 
reached  an  acceptable  level  or  when  the  receiver  de¬ 
cides  that  the  received  signal  is  of  no  interest.  Likewise, 
in  applications  where  the  receiver  is  more  interested  in 
specific  parts  of  the  signal  rather  than  the  entire  signal 
(e.g.,  content-based  image  transmission),  a  valid  ques¬ 
tion  is  how  to  transmit  and  reconstruct  a  signal  with 
different  levels  of  quality  (distortion)  in  different  tem¬ 
poral  or  spatial  regions.  For  example,  after  receiving  a 

*  Prepared  in  part  through  collaborative  participation  in 
the  Advanced  Telecommunications/Information  Distribution  Re¬ 
search  Program  ( ATIRP)  Consortium  sponsored  by  the  U.S.  Army 
Research  Laboratory  under  Cooperative  Agreement  DAAL01-96- 
2-0002  and  in  part  through  support  received  from  the  National 
Institute  of  Standards  and  Technology  grant  70NANBSH0051. 


rough  reproduction  of  a  medical  image  in  a  telemedicine 
application,  the  radiologist,  at  the  receiver,  may  want 
to  highlight  a  part  of  the  image  and  request  a  higher 
fidelity  (or  even  lossless)  replica  of  only  the  highlighted 
area. 

The  use  of  subband  decomposition  or  discrete  wavelet 
transform  (DWT)  in  image  coding  systems  has  received 
much  attention  in  recent  years  [l]-[3].  Not  only  do  these 
coding  techniques  provide  good  compression  results,  but 
also  they  are  inherently  multi-resolution.  Recently,  sev¬ 
eral  rate-scalable  subband  image  coding  systems  have 
been  proposed  in  the  literature  which  provide  very  good 
performances  [4,  5,  6,  7]. 

One  of  the  problems  in  a  wavelet-based  progressive 
transmission  scheme  is  that  the  decoder,  upon  receiving 
new  bits,  has  to  perform  the  inverse  filtering  operation 
to  reconstruct  the  image.  The  normal  approach  for  do¬ 
ing  this  is  to  apply  the  inverse  filters  on  the  decoded  ver¬ 
sions  of  all  wavelet  coefficients  to  reconstruct  a  replica 
of  the  image.  In  this  approach,  even  if  one  wavelet  coef¬ 
ficient  is  refined,  the  reconstruction  complexity  will  be 
the  same  as  in  the  case  when  all  coefficients  are  refined. 

In  a  progressive  transmission  scenario  however,  at 
each  step  of  progression,  only  some  of  the  wavelet  coef¬ 
ficients  are  refined  and  the  other  coefficients  remain  un¬ 
changed.  Therefore,  only  a  portion  of  the  image  pixels 
will  change  after  receiving  the  new  bits.  To  reduce  the 
complexity  of  reconstruction,  one  can  recompute  only 
those  pixels  of  the  image  that  need  to  be  changed.  Not 
only  does  our  proposed  scheme  provide  a  fast  recon¬ 
struction  of  the  output,  but  also  it  can  update  the  out¬ 
put  as  each  wavelet  coefficient  arrives  and  does  not  need 
to  wait  for  receiving  all  packets  before  it  starts  the  pro¬ 
cess  of  reconstruction.  This  can  be  done  without  any 
increase  in  complexity.  This  feature  of  our  reconstruc¬ 
tion  scheme  is  not  limited  to  progressive  transmission 
systems  and  can  be  used  for  any  packetized  bit  stream 
for  on-line  reconstruction  of  the  output. 
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In  Section  2,  we  propose  a  new  way  of  reconstruct¬ 
ing  the  signal  which  allows  for  updating  only  the  neces¬ 
sary  portions  of  the  signal,  instead  of  reconstructing  the 
whole  thing.  The  complexity  of  this  new  reconstruction 
approach  is  proportional  to  the  number  of  refined  coef¬ 
ficients.  Section  3  generalizes  the  approach  of  Section 
2  to  a  general  filter  bank  structure  for  subband  recon¬ 
struction  and  Section  4  contains  concluding  remarks. 

2. ONE  LEVEL  OF  DECOMPOSITION 

In  signal  coding  based  on  the  DWT,  the  input  signal  is 
typically  decomposed  into  two  components:  (i)  a  low- 
resolution  approximation  and  (ii)  a  detail  signal.  This 
results  in  decomposing  the  input  signal  into  low-pass 
and  high-pass  versions,  generally  referred  to  as  sub¬ 
bands.  Each  of  the  resulting  subbands  can  be  further 
decomposed  using  the  same  approach.  In  this  manner, 
the  DWT  decomposes  a  given  input  signal  into  a  num¬ 
ber  of  frequency  bands  [2],  At  the  receiver,  the  signal 
can  be  reconstructed  using  appropriate  inverse  DWT  fil¬ 
ters.  Fig.  1  shows  one  level  of  two-band  decomposition 
and  reconstruction  using  linear-phase,  finite-impulse  re¬ 
sponse  filter  banks  [3]. 

In  this  section,  to  convey  the  basic  idea  behind  the 
proposed  fast  reconstruction  approach,  we  limit  our¬ 
selves  to  a  one- level  decomposition  as  in  Fig.  1  where 
the  reconstructed  output  signal  can  be  written  as: 


y(n)  =  E lh=-Nh  x'h(n  ~  m)h(m)  + 


EZ9=-No  x'g(n  ~  m)g(m), 


(1) 


where 


4(n)  = 


x'Jn)  = 


xhV- 


9 ,  n  even  ,  .  .  „  ,  2  \ 

O  '2  nodd 


•r‘,(f)  neven  x'(A  =  X(72\ 
0  n  odd  ’  h 


(2) 


and  2Nh  +  1  and  2 Ng  +  1  are  the  lengths  of  the  h  and 
g  filters,  respectively. 

To  simplify  the  notation,  we  define  Xh(a)  =  xg(a)  =  0 
if  a  ^  Z.  So, 


,  n , 


x'i{n)  =  Xi{-%  i  =  h,g. 


(3) 


Using  the  above  notation,  the  output  can  be  written  as: 
y(n)  =  E%L_Nxh(^)h(m)+ 


E Zl_Nxg(^)g(m). 


(4) 


Now,  let  us  consider  a  progressive  transmission  scheme 
where  yi(n ),  l  =  1,2 ,■■■  ,K.  represents  the  different  re¬ 
finements  of  the  output.  The  refinement  of  the  out¬ 
put  can  be  achieved  by  refining  the  input  coefficients 
successively  using  b\(n)  bits/sample  to  quantize  Xj(n) 
(i  =  h,g)  at  level  l.  The  values  of  b\(n)  are  de¬ 
fined  by  a  bit  allocation  algorithm  or  are  specified  in¬ 
herently  in  the  coding  system.  Note  that  b\(n)  < 
b]+ i(n)  and  it  is  possible  -  and  crucial  to  note  - 
that  we  might  have  b\(n)  =  b\+1(n),  for  some  i.  We 
use  Xi(n ,  b\),Xi(n,  b\),  ■  ■  ■ , Xi(n,  tfK)  to  represent  the  de¬ 
coded  versions  of  the  subband  coefficients  Xj.(n),  at  dif¬ 
ferent  levels  of  refinement  (x,  (n.  b) )  is  the  coarsest  rep¬ 
resentation  of  Xi(n )  using  b\(n)  bits  for  quantization 
and  xt(n.  b'K)  is  the  finest  representation  using  blK(n) 
bits  for  quantization). 

Using  the  normal  method  of  reconstruction  in  sub¬ 
band  coding  systems,  for  each  l  =  1,2 ,  •  •  -  ,  AT,  upon 
obtaining  x,fL  and  xfj .  the  decoder  reconstructs  the  l^1 
refinement  of  the  output,  yi ,  according  to: 


yi(n)  =  E1ZL- k  xh(^M)h(m)+ 


^  ll h _ 

2^m=-Ng  Xg\  2  5  ul 


2  ’  “ l 

b9i)g{m) 


(5) 


In  other  words,  given  that  yi  is  reconstructed  at  the 
decoder,  upon  receiving  additional  bits  from  which 
Xh(n,  6f+1)  and  xg(n,b9l+1)  -  the  (l  +  l)st  refinement  of 
x fr  and  xg  -  are  obtained,  the  inverse  DWT  filters  are 
used  anew  to  compute  an  updated  replica  of  the  signal 
-  the  ( l  +  l)st  refinement  of  y. 

In  this  paper,  we  propose  an  alternative  method 
in  which  we  define  A yi(n)  =  yi+i(n)  —  yi(n)  and 
A Xi{n,b\)  =  Xi(n,  b\+l)  -  Xi{n,b\),  i  =  h,g.  Clearly, 
A  yi,  l  =  1,  2f>  ■  •  ,  K ,  can  be  computed  as  follows: 


&yi{n)  =  EnL-N,  A xh{^p-M)h(m)  + 


,b9i)g{m). 


(6) 


Then,  to  reconstruct  a  new  refinement  of  the  signal, 
yi- |_i,  it  suffices  to  use 


Vi+  i(n)  =  yi{n)  +  A  y/(n). 


(7) 


Now  let  us  consider  the  complexity  implication  of  this 
new  approach  and  its  comparison  with  the  normal  ap¬ 
proach  using  (5).  If  the  number  of  samples  of  the  input 
signal  (x)  is  N,  then  the  number  of  samples  in  x/,  and 
xg  will  be  y .  If,  at  each  level  of  refinement  l,  we  use  (5) 
to  compute  yi,  it  costs  L  =  -y(2A Th  +  1)  +  -j(2 Ng  +  1) 
multiplications,  a  quantity  which  is  independent  of  the 
number  of  refined  samples.  Therefore,  after  K  levels  of 
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refinement  the  computational  complexity  is  KL  multi¬ 
plications.  Using  (6)  and  (7)  to  update  the  output  signal 
costs  2 Nj  + 1  multiplications  per  each  nonzero  Ax,  coef¬ 
ficient,  i  =  h,g  (number  of  nonzero  Ax,  <  -y)-  If  °il  y 
wavelet  coefficients  [x,,  i  =  h,g)  are  refined  at  the  l^1 
level  of  refinement,  then  the  total  number  of  multipli¬ 
cations  using  (6)  and  (7)  to  update  the  output  signal  is 
Yi=  l  t^hl (2iV/j  +  1)  +  Yi= l  aglY^Ng  +  1).  Note  that 
for  large  K' s,  Yi= i  ail  <<  K  f°r  existing  rate-scalable 
image  coding  systems  [4,  5,  6,  7]  and  therefore  the  com¬ 
plexity  of  the  proposed  approach  is  much  less  than  the 
complexity  of  the  normal  approach.  For  example,  in  the 
embedded  zerotrees  wavelet  (EZW)  coding  approach  of 

[4] ,  not  only  are  there  many  insignificant  coefficients  at 
low  bit  rates  which  are  not  transmitted  and  are  assumed 
to  be  zero,  but  also  the  significant  coefficients  are  refined 
in  dominant  and  subordinate  passes  (at  each  pass  some 
of  the  “significant”  coefficients  remain  unchanged). 

In  addition  to  reducing  the  complexity  of  reconstruc¬ 
tion,  using  (6)  and  (7)  to  reconstruct  the  output  pro¬ 
vides  the  capability  to  update  the  output  signal  upon 
receiving  the  new  refinement  of  each  wavelet  coefficient, 
thus  allowing  for  an  on-line  update  of  the  output.  Typi¬ 
cally,  for  any  signal  coding  system  the  output  bit  stream 
is  transmitted  in  the  form  of  packets.  Consider  a  gen¬ 
eral  signal  coding  system  (or  one  level  of  refinement  of 
a  progressive  system)  and  assume  that  the  output  is 
transmitted  using  M  packets.  The  size  of  each  packet 
and  thus  the  value  of  M  depends  on  the  network  char¬ 
acteristics.  At  the  receiver,  either  we  should  wait  to 
receive  all  M  packets  or  we  can  reconstruct  M  inter¬ 
mediate  versions  of  the  output  -  one  for  each  newly  re¬ 
ceived  packet.  Using  the  normal  reconstruction  formula 

(5) ,  results  in  a  delay  in  the  former  case  and  an  increase 
in  computational  complexity  by  a  factor  of  M  in  the 
latter  case.  Using  (6)  and  (7)  to  reconstruct  the  out¬ 
put  will  provide  M  intermediate  versions  of  the  output 
without  increasing  the  complexity  or  delay.  Also,  note 
that  if  after  receiving  a  rough  reproduction  of  the  signal, 
we  just  need  to  refine  a  portion  of  wavelet  coefficients 
(corresponding  to  a  spatial  location),  the  computational 
cost  using  (6)  and  (7)  is  proportional  to  the  number  of 
transmitted  coefficients  (filter  inputs)  and  less  than  the 
computational  cost  of  the  normal  scheme. 

To  update  the  output,  we  need  to  know  the  set  of 
output  samples  whose  values  are  influenced  by  the  re¬ 
finement  of  each  coefficient.  This  set  depends  on  the 
location  of  the  refined  coefficient.  Also,  the  samples 
near  the  boundaries  need  a  special  treatment  because  of 


the  symmetric  extension  of  the  input  signal1.  Receiv¬ 
ing  Axh(k)  for  the  coefficients  sufficiently  far  from  the 
boundaries,  -y-  <  k  <  ;  affects  only  the  output 

values  y{2k  +  n),  n  =  — N —  (TV/,,  —  1),  •  •  •  ,  Nf,.  When 
the  received  coefficient  is  near  the  origin,  0  <  k  < 
the  output  samples  yin),  n  =  0, 1,  •  •  ■  ,  Nf,  +  2k,  are  af¬ 
fected.  For  the  remaining  coefficients,  A  ~^h  <  k  <  - y, 
the  output  values  y{2k  +  n),  n  =  —N^,  ■  ■  ■  ,  N  —  2k,  are 
affected.  Replacing  Xf,  with  xg  and  Nf,  with  Ng  provides 
the  same  indices  for  the  filter  g  and  coefficients  xg. 

The  same  procedure  can  be  applied  to  more  than  one- 
level  of  decomposed  signals.  The  two-dimensional  (2-D) 
extension  of  the  approach  is  also  straightforward  when 
we  use  a  separable  2-D  transform. 

3. GENERAL  FILTER  BANK 

In  this  section,  we  extend  the  idea  for  complexity  re¬ 
duction  developed  in  the  previous  section  to  general  fil¬ 
ter  banks.  Fig.  2  illustrates  a  general  decomposition 
and  reconstruction  structure  using  a  set  of  passband 
filters  called  filter  banks.  This  structure  can  be  used 
to  perform  linear  transformations  like  the  discrete  co¬ 
sine  transform,  lapped  orthogonal  transform,  Laplacian 
pyramid,  Gabor  transform,  quadrature  mirror  filters, 
and  DWT  [8].  Although  there  exist  “fast  implemen¬ 
tations”  for  many  of  these  transformations,  the  study 
of  this  general  structure  is  worthy  of  consideration  as 
it  unifies  the  application  of  our  approach  for  different 
transformations  and  provides  a  general  methodology  for 
progressive  transmission  situations.  Furthermore,  our 
approach  does  not  depend  on  the  specific  transforma¬ 
tion  used.  In  fact,  the  ratio  of  the  computational  com¬ 
plexity  of  the  proposed  approach  to  that  of  the  normal 
approach  for  reconstruction  is  fixed  for  any  transforma¬ 
tion.  So,  the  choice  of  the  transformation  affects  the 
complexity  of  both  approaches  in  the  same  manner. 

To  make  the  argument  more  precise,  first,  let  us  con¬ 
sider  the  example  illustrated  in  Fig.  3  consisting  of  three 
levels  of  the  two-band  structure  in  Fig.  1.  Later,  we 
generalize  our  fast  reconstruction  approach  to  the  filter 
bank  shown  in  Fig.  2.  The  example  in  Fig.  3  is  equiv¬ 
alent  to  Fig.  2  when  k\  =  2,  k->  =  4,  kz  =  Aq  =  8  and 
filters  gi  and  hi,  i  =  1,2,  3, 4,  are  chosen  appropriately. 
Given  filters  h  and  g  in  Fig.  3,  the  corresponding  recon¬ 
struction  filters  in  Fig.  2  can  be  found.  The  main  equa¬ 
tion  governing  the  relationship  between  gi,  i  =  1,  2,  3,4 
and  g  and  h  is  the  equivalence  of  the  two  structures  in 

lrTo  have  a  smooth  boundary  after  reconstruction,  the  signal  is 
symmetrically  extended  near  the  boundaries. 
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Fig.  4  [8].  Using  Fig.  4,  one  can  establish  the  following 
equations  relating  different  filters  in  Figs.  2  and  3: 


=  G(^), 

(8) 

g2(*) 

=  G(^)F7(^), 

(9) 

g3(^) 

=  G(z4)H(z2)H(z)1 

(10) 

g4(^) 

=  H(z4)H(z2)H(z). 

(ID 

Although  the  reconstruction  operation  in  Figs.  2  and 
3  are  the  same,  the  number  of  multiplications  in  Fig.  2 
is  more  than  the  number  of  multiplications  in  Fig.  3. 
Specifically  the  output  signal  in  Fig.  2  can  be  written 
as: 


N  i 

y{n)  =  Y  x'^n  ~  + 

m=  —  Nj 
N2 

Y  x'2(n  -  m)g2(m)  + 

m=  —  N2 

n3 

Y  x'3[n-in)gz{in)  + 

m=—N3 

N.i 

Y  x'4  (n-m)g4(m),  (12) 

TO  —  —  Ni 


where 

x[{n) 

x'2(n)  = 

x'z  (n)  = 

x4  (n)  = 


x\(^)  n  even 
0  n  odd 


®2( 

-1 

4  1 

II 

e 

0 

otherwise 

3^3  ( 

t) 

00 

II 

e 

0 

otherwise 

a:4 

(f) 

00 

II 

e 

0 

otherwise 

X[(z)=X1(z2), 
X'(z)=X2(z4), 
X'(z)=X3(z% 
X'(z)=X4(z%  (13) 


and  2Nt  +  1  is  the  length  of  gj  and  Ni  =  Ng,  N2  = 
2 Ng  +  Nh,  N3  =  ANg  +  3 Nh,  and  N4  =  7Nh.  So,  using 
all  coefficients  to  reconstruct  the  output,  we  need  L\  = 
N\?XLk±  +  2Mtl  +  M|±l  +  Mi+1]  multiplications  in 
Fig.  2  and  L2  =  N[(2Nh  + 1)(|  +  \  +  I)  +  (2 Ng  + 1)(±  + 
\  +  ^)]  multiplications  in  Fig.  3.  After  computing  L\  in 
terms  of  AT/,,  and  Ng ,  we  can  show  that 

Li  =  N[l+3Ng+3Nh]  >L2  =  N[^(Nh+Ng  + 1)],  (14) 

establishing  that  the  number  of  multiplications  corre¬ 
sponding  to  the  implementation  of  Fig.  3  is  less  than 
that  of  Fig.  2. 


Now,  we  concentrate  on  the  main  issue  of  this  section 
and  show  that  an  approach  similar  to  what  was  pro¬ 
posed  in  Section  2  can  be  used  for  the  filter  bank  shown 
in  Fig.  2.  To  simplify  the  notation,  we  define  Xj\a)  =  0 
\i  a  tfiTL.  So, 

Tl  Tl 

x'i(n)  =  Xi(  — ),  *  =  1,2,3,  and  x4(n)  =  ®4(-).  (15) 

Using  the  above  notation,  the  output  can  be  written  as: 

Ni 

,  x  v — >  /  Tl  771 .  .  . 

y(n)=  Y  Xl^  2 — )^i(m)  + 

m=—Ni 

n2 

E.  Tl  171 .  .  . 

X2(  4  )g2(m)  + 

m—  x2 

n3 

vA  .n- in.  ,  . 

Y  x^~ — >93(m)  + 
m=—N3 

Y  -4  (^pW).  (16) 

m=—N4 

In  the  normal  approach  for  reconstructing  the  output 
of  a  progressive  transmission  scheme,  the  different  lev¬ 
els  of  refinement  of  the  output,  y4,  y2l  ■  ■  • ,  y k ,  are  com¬ 
puted  after  receiving  Xi(n,  b\),Xi(n,  b'2).  •  •  •  . .r,;(n.  b‘K)  as 
follows: 

Ai  _ 

yi(n)=  Y  Xl ( n  2  m  ’  b^gi + 

m——Ni  “ 

N-2 

E.  Tl  171  9\  /  \ 

(  4  ,b,)g2(in)  + 

m  A", 

N3 

vA  .  n  -  in  ,  o .  .  . 

X,  g3(  8  ,ftf)g3(m)  + 

m=—N3 

N4  _ 

Y  M)gA{m).  (17) 

TO=  — 

As  before,  we  define  A yi(n)  =  yi+\(n)  —  yi{n)  and 
A Xi{n,b\)  =  Xi(n,  b\+l)  -  Xi{n,b\),  where  A yu  l  = 
1,2,  ■  -’,K,  is  given  by: 

Ni  _ 

Ayi  (n)  =  X  AaA  ( n  2?n )  )9i  (m)  + 

m=—Ni 

n2 

A  ,n  —  m  o  \  /  x 

X  Ax2(  ^  ,bi)g2(m)  + 

m——N2 

N3 

a  .n  —  m  ,o,  .  . 

2Z  Aa;3(  8  ,5f)g3(m)  + 

III  :V;> 

*V4  _ 

X  Ax4  ( n  m,bf)g4{in).  (18) 

•/)/  :V.| 
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Then,  given  that  yi(n)  is  available  at  the  decoder, 
yi+i(n)  can  be  computed  by  adding  A yi(n)  to  yi{n). 
Equation  (18)  is  a  generalization  of  (6)  and  therefore 
provides  a  fast  reconstruction  scheme  as  was  discussed 
in  Section  2. 

To  update  the  output,  we  need  to  know  the  list  of 
output  samples  which  must  be  updated  due  to  receiving 
each  coefficient.  Table  1  shows  the  required  updates  for 
the  coefficients  that  are  sufficiently  far  from  the  bound¬ 
aries.  The  first  row  contains  the  filter  coefficients  which 
are  used  to  update  the  corresponding  outputs  in  the 
second  row.  To  make  the  details  of  the  reconstruction 
scheme  more  clear,  let  us  provide  a  numerical  example. 

3.1. EXAMPLE 

Although  our  scheme  works  for  any  scalable  transform 
coding  system,  we  pick  a  particular  coding  system  to 
show  the  practical  usefulness  of  our  scheme.  We  use  the 
EZW  approach  of  [4]  for  this  purpose.  In  [4],  an  8  x  8 
wavelet  decomposed  image  has  been  used  to  present  the 
details  of  the  EZW  algorithm.  We  use  the  same  exam¬ 
ple  (Table  2)  and  assume  that  refining  each  coefficient 
affects  all  output  pixels  in  the  same  row  or  column  (a 
logical  assumption  on  the  length  of  filters  for  an  8  x  8 
image).  The  normal  reconstruction  approach  takes  1024 
multiplications  to  create  each  refinement  of  the  output. 
Without  going  through  the  details  of  the  EZW,  as  it 
is  shown  in  [4],  the  first  dominant  pass  results  in  re¬ 
fining  coefficients  in  locations  11,  12,  13,  and  54  (the 
numbers  indicate  the  row  and  column  indices,  respec¬ 
tively).  Reconstructing  the  image  using  our  approach 
takes  160  multiplications  instead  of  1024.  The  first  sub¬ 
ordinate  pass  (the  second  level  of  refinement)  refines  the 
same  four  coefficients  [4]  and  thus  results  in  the  same 
amount  of  computational  reduction.  The  second  dom¬ 
inant  pass  (the  third  level  of  refinement)  changes  coef¬ 
ficients  in  locations  21  and  22.  Our  approach  requires 
80  multiplications  to  reconstruct  the  third  image.  The 
second  subordinate  pass  (the  fourth  level  of  refinement) 
refines  the  aforementioned  six  wavelet  coefficients.  To 
create  the  fourth  image,  our  approach  requires  240  mul¬ 
tiplications.  Therefore,  to  reconstruct  the  first  four  im¬ 
ages,  the  fast  scheme  needs  a  total  of  640  multiplications 
compared  to  4096  multiplications.  Also,  note  that  since 
we  have  had  16  coefficient  refinements  in  the  first  four 
passes,  we  could  have  reconstructed  16  intermediate  im¬ 
ages  without  any  additional  computational  complexity. 
Constructing  these  16  images  would  cost  16384  multipli¬ 
cations  using  normal  wavelet  reconstruction  formulas. 

4. CONCLUSIONS 


We  have  proposed  a  scheme  for  fast  reconstruction  of  a 
subband-decomposed  progressively  transmitted  signal. 
Using  the  proposed  approach,  we  can  update  the  image 
after  receiving  each  new  coefficient  and  create  a  contin¬ 
uously  refined  perception  without  any  extra  computa¬ 
tional  cost  (compared  to  the  case  that  we  reconstruct 
the  image  after  receiving  each  packet).  In  existing  scal¬ 
able  image  coding  systems  [4]- [7],  insignificant  coeffi¬ 
cients  remain  zero  at  the  decoder  for  a  few  steps  of  re¬ 
finement.  Also,  at  each  step  of  the  refinement,  some 
of  the  coefficients  remain  unchanged.  Unrefined  coeffi¬ 
cients  do  not  add  to  the  computational  complexity  of 
our  new  approach  for  reconstructing  the  image.  On  the 
other  hand,  the  complexity  of  the  normal  reconstruction 
scheme  does  not  depend  on  the  value  of  the  coefficients 
(or  refinements). 

In  some  applications,  after  receiving  a  preliminary 
draft  of  the  image,  the  receiver  only  needs  the  upgrade 
of  a  particular  portion  of  the  image.  Using  our  proposed 
approach,  the  complexity  reduction  in  this  case  is  equal 
to  the  ratio  of  the  number  of  pixels  in  the  required  area 
to  the  number  of  pixels  in  the  whole  image. 

Future  work  includes  the  generalization  of  the  scheme 
for  filter  banks  with  integer  coefficients  which  replace 
multiplications  by  binary  shifts.  2 
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Figure  1:  Two-Band  Subband  Analysis  and  Synthesis 
Structure. 
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Figure  2:  Subband  Analysis  and  Synthesis  Filter  Banks. 


x\(k) 

gV-Ni) 

■■■  5i (0) 

5i(JVi) 

NI  <k<  N~0Nl 

y(2k-Nl) 

•  •  •  1/(2*) 

•••  //(2  A-  +  .V, ) 

X-2  (k) 

52  (— iV2) 

•••  52(0) 

52  (Ah) 

<  k  <  N-,^ 

y{4k-N2) 

•  •  •  2/(4*) 

//(I  A-  +  \%) 

X3  (k) 

53  (—Ah) 

•••  5s(0) 

53(W) 

-V.  <  /,  . 

y(8k-N3) 

•  •  •  2/(8*) 

5(8*  +  Ah) 

xa  {k) 

54  (-Ah) 

•••  54(0) 

54  (W) 

N*  <  k  <  N~XN 4 

y(8k  -  Ah) 

•  •  •  2/(8*) 

•••  y(8k  +  Ni) 

Table  1:  Update  Table  for  Coefficients  Sufficiently  Far 
from  the  Boundaries. 
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-1 
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3 

-12 
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-7 
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9 

-9 

-7 

-14 

8 

4 

-2 

3 

2 

-5 

9 

-1 

47 

4 

6 

-2 

~2 

3 

0 

-3 

2 

3 

-2 

0 

4 

2 

-3 

6 

-4 

3 

6 

3 

6 

5 

11 

5 

6 

0 

3 

-4 

4 

Table  2:  Example  of  3-level  wavelet  transform  of  an 
image  [4], 
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